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Improving  foF2  Prediction  for  the  Sunrise 

Transition  period 


I.  IM'HOIH  CTION 

During  the  transition  period  of  sunrise  (and  sunse*),  the  critical  frequency 
foF2  of  the  ionospheric  F*layer  changes  at  a  rapid  rate,  demanding  more  frequent 
changes  in  the  operating  frequency  of  an  OTH-B  radar  to  maintain  a  relatively 
constant  barrier  location.  The  observations  as  well  as  predictions  from  iono¬ 
spheric  models  show  that  the  change  in  foF2  is  between  20  and  60  percent  (average 
3C  percent)  per  hour  during  the  sunrise  period  (05-09  LT)  from  October  through 
March  independent  of  the  phase  of  the  solar  activity  cycle. 

An  OTH-B  surveillance  radar  is  required  to  maintain  a  barrier  at  a  distance 
of  20C0km  from  the  radar.  During  regular  daytime  and  nighttime  operations, 
when  the  ionosphere  is  in  general  a  horizontally  stratified  medium,  a  backscatter 
system  can  Illuminate  a  large  area  beyond  the  skip  distance.  During  sunrise  hours, 
however,  the  sunrise-related  tilt  of  the  F-layer  results  in  a  narrowing  of  the  back¬ 
scatter  radar  "footprint,"  that  is,  the  useful  range  for  target  detection.  In  this 
discussion,  we  consider  a  nominal  value  of  500  km  (275  nm)  as  the  morning  hour 
useful  barrier  width.  In  order  to  track  targets  crossing  this  narrow  barrier,  it  is 
considered  important  that  the  radar  footprints  before  and  after  an  operational  fre¬ 
quency  change  overlap  by  approximately  50  percent  or  250  km.  This  will  ensure 
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that  about  50  percent  of  the  already  established  targets  will  be  Illuminated  after 
the  frequency  switch,  providing  track  continuity  important  to  the  identification  pro¬ 
cess  during  times  when  aircraft  on  the  average  cross  the  barrier  in  approximately 
35  minutes. 

Therefore,  during  tlie  sunrise  period,  the  frequency  would  have  to  be  changed 
every  time  the  bander  moves  towards  the  radar  by  250km.  At  a  distance  of 
2250km  from  the  transmitter,  that  is,  in  the  center  of  a  nominal  barrier  extending 
from  2C00km  to  25C0km,  the  M-factor  (which  relates  the  radar  operation  fre¬ 
quency  to  foF2  at  the  reflection  point)  changes  at  a  rate  of  2.  3  percent  for  every 
100  km.  Conversely,  this  means  that,  for  every  2,3  percent  change  in  M-factor 
(and  therefore  in  foF2),  the  barrier  will  move  100  km  towards  the  radar  during 
sunrise.  Therefore,  for  every  5.8  percent  change  in  M-factor  or  foF2,  the  bar¬ 
rier  will  shift  the  maximum  permissible  250  km.  To  keep  up  with  a  30  percent 
per  hour  change  in  foF2  during  sunrise  hours,  and  maintain  a  50  percent  overlap 
of  the  barrier  through  each  frequency  change,  the  operating  frequency  of  the  radar 
has  to  be  adjusted  at  least  5  times  per  hour.  It  is  therefore  desirable  to  have  re¬ 
liable  short  term  predictions  of  foF2  for  the  coverage  region  of  interest,  to  facili¬ 
tate  t!;e  necessary  rapid  changes  in  operating  frequency. 

Inspection  of  throe  readily  available  techniques  showed  that  they  would  not 
solve  the  problem.  These  techniques  are:  the  use  of  existing  models,  the  extrap¬ 
olation  of  foF2  from  measurements  at  a  different  location  (such  as  the  radar  site), 

and  extrapolation  of  the  expected  sunrise  transition  from  the  foF2  sunrise  transi- 
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tion  measured  earlier  at  a  more  easterly  location.  The  lTS-78  or  lONCAP 

2 

monthly  median  prediction  models  (Lloyd  et  al  )  yield  an  foF2  prediction  error  of 
±1.  5  MHz,  These  models  are  unable  to  provide  the  required  accuracy  of  ^  0, 5  MHz 
needed  by  the  real  time  operational  systems  (this  requirement  for  sunrise  transi- 

3 

tion  will  be  substantiated  in  Section  2).  Rush  and  Miller  have  shown  that  the  two 
station  extrapolation  scheme  for  foF2  prediction  breaks  down  (correlation  coeffi¬ 
cient  <  0.7)  beyond  1800  km  in  the  E-W  direction  and  beyond  1000  km  in  the  N-S 
direction.  The  distance  for  which  the  correlation  coefficient  equals  0.7  is  subse- 


1.  Barghausen,  A.L.,  Finney,  J.  W. ,  Proctor,  L.  L. ,  and  Schultz,  L.D.  (1969) 
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2.  Lloyd,  J.L.,  Haydon,  G.W.,  Lucas,  D.L.,  and  Teters,  L.R.  (1978)  Estimat¬ 
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qucntly  referred  to  as  the  correlation  distance.  Davis  and  Rush  studied  correla¬ 
tions  using  European  and  American  stations  as  well  as  the  net  of  Russian  iono¬ 
spheric  stations  for  prediction  of  foF2  during  the  sunrise  transition  period.  They 
conclude  that  the  extrapolation  scheme  of  predicting  foF2  for  the  OTH-B  coverage 
area  on  the  East  Coast  of  the  United  States  from  observations  at  European  stations 
does  not  work.  There  is  an  additional  complication  because  the  reflection  points 
of  interest  to  the  OTH-B  coverage  and  the  corresponding  European  stations  are  at 
the  same  geographic  (solar  Influence)  latitude  but  are  at  quite  different  geomagnetic 
(auroral/trough  influence)  latitudes. 

The  ionospheric  reflection  region  of  an  OTH-B  radar  is  typically  10C0-1250km 
away  from  the  transmitter.  For  spatial  extrapolation  from  the  radar  ionosonde  to 
tiie  ionospheric  reflection  region,  tiiis  separation  distance  Ijarely  satisfies  the  cor- 

3 

relation  range  determined  by  Rush  and  Miller  depending  on  the  direction  (N-S, 
E-W)  of  radar  illumination.  In  addition,  the  sunrise  observations  at  the  East 
Coast  Radar  Site  v/ill  be  at  the  same  time  (along  the  N-S  meridian),  or  at  a  later 
time,  but  never  earlier  than  sunrise  in  the  OTH-B  reflection  regions.  Thus,  even 
for  real  time  specification,  the  vertical  ionosonde  observations  from  the  transmit¬ 
ter  site  have  to  be  extrapolated  both  in  time  as  well  as  in  distance  in  applying  the 
two  station  method  for  OTH-B  coverage.  Prediction  must  still  be  made  for  the 
next  hour(s).  The  result  of  such  a  scheme  may  be  acceptable  for  non-transition 

periods  but  is  not  useful  for  the  transition  period  in  the  view  of  the  Rush  and 

3  4 

Miller  and  Davis  and  Rush  investigations. 

To  obtain  foF2  measurements  for  better  real-time  specification  of  the  iono¬ 
sphere,  and  other  requirements  such  as  the  reliable  h'(f)  determination  at  key  lo¬ 
cations,  ionosondes  have  been  deployed  at  Goose  Bay  (GS),  Labrador,  and  Argen- 
tia  (AG),  Newfoundland,  and  deployment  is  planned  for  ionosondes  at  Bermuda  (BR) 
and  Wallops  Island,  Va  (WA).  These  stations  will  provide  real-time  ionospheric 
data  from  under  or  close  to  the  OTH  midpoint  areas  for  tlie  East  Coast  system. 

(Sec  locations  with  respect  to  the  midpoint  ionosphere  shown  in  Figure  1. )  This 
will  eliminate  or  greatly  reduce  the  need  for  spatial  extrapolation  of  the  observa¬ 
tions,  making  these  data  very  reliable  for  updating  of  the  OTH-B  virtual  height/ - 
coordinate  conversioti  tables.  However,  for  the  period  of  fast  foF2  changes,  a 
better  short  term  prediction  scheme  is  still  needed  to  guide  the  frequency  manager 
in  the  necessary  rapid  frequency  selection.  Therefore,  this  study  was  conducted 
to  find  a  suitable  approach  of  predicting  foF2  from  prior  days'  and  real  time  foF2 


4.  Davis,  R.M,,  Jr.,  and  Rush,  C.  M.  (19S3)  Feasibility  of  Forecasting  foF2 

Disturbances  During  the  Sunrise  Transition  Period,  Technical  Memorandum 
Series  NTIA -TM-33-87,  U.S,  Department  of  Commerce,  National  Telecom¬ 
munications  and  Information  Administration. 
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Figure  1.  OTH  Coverage  Area  of  the  East  and  West  Coast  Based  Radar  Systems 


observations  at  these  stations  within  the  OTH  midpoint  area.  This  report  presents 
the  results  of  the  study. 

Section  2  defines  in  detail  the  problem  encountered  in  radar  frequency  manage¬ 
ment  during  the  sunrise  period.  Section  3  presents  the  data  base  used  in  this 
study.  The  different  studies,  such  as  the  accuracy  of  lONCAP  predictions,  the 

4 

tv/o  station  method  of  Davis  and  Rush,  and  the  study  of  the  scheme  proposed  in 
this  paper  did  not  always  use  the  same  data  base.  In  Table  1,  we  identify  the  sta¬ 
tions  used  for  the  present  study  as  well  as  those  used,  in  addition  or  instead,  in 
the  other  studies. 

The  analysis  in  Section  4  presents  several  new  approaches  to  short  term  foF2 
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Table  1.  Ionospheric  Stations  Used  in  the  Analysis 
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predictions  that  were  evaluated  for  usefulness,  such  as  the  use  of  prior  clays'  ave¬ 
rage  foF2  at  the  corresponding  LT,  or  the  combination  of  previous' day(s)  hourly 
slope  of  foF2  with  the  last  ol  served  foFS  for  a  1-hr  prediction.  As  a  baseline,  to 
judge  improvement  of  the  short  term  predictions,  we  determined  prediction  errors 

O 

using  the  lONCAF’  model'  alone. 

This  analysis  shows  that  the  best  prediction  for  a  1-hr  interval  is  the  combi¬ 
nation  of  '.he  currently  observed  foF2,  the  average  hourly  slope  determined  using 
foF2  measuremcr.ts  for  the  past  four  days,  and  a  correction  in  the  slope  as  a  I'unc- 
'ion  of  the  deviation  of  the  oi;sorved  fcjF2  from  the  past  4-day  average  foF2  value 
for  the  particular  hour.  The  procedure  for  such  prediction  of  foF2  is  summarized 
in  t!ie  conclusion. 

The  Experimental  Radar  System  (ERS)  consistently  experienced  frequency 
management  difficulties  during  2-3  hr  of  the  sunrise  period,  v/ith  the  system  per¬ 
forming  poorly  for  8  to  12  percent  of  the  operation  time  on  a  routine  basis.  It  is 
essential  to  realize  that  these  observed  difficulties  occurring  8  to  12  percent  of  the 
time  were  not  a  random  problem.  The  improvements  possible  with  the  proposed 
prediction  technique  should  increase  the  overall  performance  of  the  system  and  at 
least  randomize  poor  performance  in  the  difficult  sunrise  period. 
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2.  DKFIMTION  OF  THK  SI  NHISK  foF2  PKOBLEM 


The  foF2  contours  for  several  stations  of  OTH-B  interest  are  shown  in  the  fol¬ 
lowing  figures  covering  high  and  low  phases  of  sunspot  activity.  The  stations  se¬ 
lected  are  Goose  Bay,  St.  Johns,  Ottawa,  Wallops,  and  Argucllo.  In  the  OTH-B 
coverage  area  shown  in  Figure  1,  Goose  Bay  is  centered  in  the  northern  3C®  beam 
sector  of  the  East  Coast  Radar  System  (ECRS).  It  is  located  at  the  equatorward 
boundary  of  the  auroral  oval.  The  station  moves  in  and  out  of  the  auroral  region 
due  to  rotation  of  Earth.  St.  Johns  is  south  of  the  auroral  zone  and,  at  night,  lies 
close  to  or  within  the  mid-latitude  F-layer  trough  region.  Ottawa  is  at  the  same 
C.G.  latitude  as  St.  Johns,  lies  west  of  the  radar  coverage,  is  used  to  confirm  the 
ii>nospheric  behavior  at  St.  Johns,  and  replaces  St.  Johns  whenever  St.  Johns  data 
are  lacking.  O  tavva  is  also  located  at  a  C.G.  latitude  similar  to  Sitka,  Alaska,  a 
sounder  site  selected  for  the  West  Coast  Radar  Systems  (WCRS).  Wallops  repre¬ 
sents  the  center  of  the  southern  coverage  area  of  180°  from  the  ECRS  and  is  unaf¬ 
fected  by  auroral  or-  ionospheric  trough  phenomena.  Arguello  represents  the  iono¬ 
spheric  conditions  for  large  parts  of  the  coverage  region  of  the  WCRS  and  for  the 
southernmost  beams  of  the  ECRS.  Thus,  these  stations  cover  most  of  the  mid¬ 
point  ionospheres  relevant  to  the  east  and  west  coast  radars. 

The  foF2  contours  from  the  monthly  median  observations  and  from  the 
lONCAP  monthly  predictions  are  shown  for  Goose  Bay,  St.  Johns,  Ottawa,  Wal¬ 
lops,  and  Arguello  in  Figures  2a-2e,  for  the  high  sunspot  activity  (SSN  =  160)  pe¬ 
riod  from  June  1979  to  June  1980.  Figures  3a-3d  cover  the  low  sunspot  activity 

(SSN  =  14)  period  for  the  calendar  year  1976,  No  data  are  available  from  Arguello 

2 

in  197o,  In  each  figure,  the  bold  lines  show  the  lONCAP  monthly  median  predic¬ 
tions  (using  the  observed  smoothed  sunspot  number),  and  the  dashed  lines  show  the 
monthly  median  observations.  The  monthly  median  observations  are  derived  from 
the  routine  hourly  ionospheric  soundings  from  the  stations. 

The  features  common  to  all  these  figures  are:  (1)  a  good  agreement  (v/ithin 
±1  Mliz)  between  predictions  and  observations;  (2)  the  bunching  of  the  contours  in 
winter  months  both  predicted  and  observed  in  the  06-08  LT  time  window  showing  a 
rapid  increase  in  foF2  during  the  sunrise  hours;  (3)  a  factor  of  2  increase  in  noon 
foF2  from  sunspot  minimum  to  sunspot  maximum;  and  (4)  the  presence  of  a  strong 
sunrise  transition  in  foF2  for  all  the  v.inter  months  independent  of  the  phase  of  the 
solar  cycle  activity.  The  figures  show  that  although  the  models  reproduce  most 
of  the  observed  foF2  behavior,  they  are  unable,  even  in  the  median,  to  provide  the 
prediction  accuracy  of  <  C.SMIIz  needed  in  real  time  by  the  operational  systems 
during  sunrise  periods.  The  general  accuracy  of  lONCAP  monthly  median  foF2s 
is  il  MHz,  and  somewhat  worse  for  sunrise  transition  periods,  when  compared  to 
the  median  of  the  observed  data. 
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Figure  2a.  Monthly  Median  Contours  of  foF2  for  High  Sunspot  Activity  (SSN  =  160)  for  Goose  Bay 


Figure  2b.  Monthly  Median  Contours  of  foF2  for  High  Sunspot  Activity  (SSN  =  160)  for  St.  Johns 
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Figure  2c.  Monthly  Median  Contours  of  foF2  for  High  Sunspot  Activity  (SSN  =  160)  for  Ottawa 
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Figure  2e.  Monthly  Median  Contours  of  foF2  for  High  Sunspot  Activity  (SSN  =  160)  for  Arguello 


9Z6I 
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Figure  3a.  Monthly  Median  Contours  of  foK2  for  Low  Sunspot  Activity  (SSN  *  14)  for  Goose  Bay 


9Z6I 
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Figure  3b.  Monthly  Median  Contours  of  foF2  for  Low  Sunspot  Activity  (SSN  =  14)  for  St.  Johns 


9Z6I 
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Figure  3c.  Monthly  Median  Contours  of  foF2  for  Low  Sunspot  Activity  (SGM  =  14)  for  Ottawa 
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Figure  3d.  Monthly  Median  Contours  of  foF2  for  Low  Sunspot  Acti/ity  (SSN  =  14)  for  Wallops 


For  determining  the  rapid  change  in  foF2  (from  Figures  2  and  3)  during  sun¬ 
rise,  a  quantitative  measure  Pj  is  defined  as  the  percentage  in  change  of  foF2  over 
1  hr. 


(foF2i  +  1  -  foF2i) 

P.  =  200  X  -  (1) 

1  (foF2i  ^  1  +  foF2i) 

where  foF2.  refers  to  the  hourly  observations  at  a  station  at  the  ith  hour.  This 
percent  change  in  foF2  is  shown  for  the  sunrise  period  in  Figure  4  for  high  sunspot 
activity  (SSX  =  138)  and  in  Figure  5  for  low  sunspot  activity  (SSN  =  17), 

These  figures  show  that  during  sunrise,  gradients  >  20  percent  last  for  about 
3  hr.  This  rapid  increase  in  foF2  is  present  for  at  least  8  months  (centered  on  win¬ 
ter)  of  the  year.  Gradients  as  high  as  60  percent  per  hour  for  short  periods  are 
present  in  some  months.  In  years  of  low  sunspot  activity.  Goose  Bay  (Figure  5) 
also  shows  a  30  percent  gradient  in  foF2  during  sunset  hours  in  the  winter  months. 
None  of  the  other  stations  show  a  similar  rapid  foF2  change  during  sunset. 

The  radar  operation  frequency  is  related  to  the  midpoint  foF2  by  the  equation, 

^RADAR  ^  ^ 

where  the  multiplying  factor  M  is  a  function  of  h^  and  y^  (the  height  and  thick¬ 
ness  of  the  layer)  and  the  range  to  the  barrier.  To  accommodate  the  rapid  in¬ 
crease  in  foF2  during  sunrise,  the  operating  frequency  of  the  radar  has  to  be 
changed  in  step  with  the  sunrise  changes  in  foF2.  The  factor  relevant  to  the  fre¬ 
quency  management  is  the  M  factor,  which  changes  with  skip  distance.  To  sim¬ 
plify  the  discussion,  we  will  consider  a  parabolic  layer  shape,  an  assumption  used 
in  the  computations  presented  here  (Appleton  and  Beynon^).  While  other  layer 
shapes  may  slightly  alter  the  M  factor  for  a  given  skip  distance,  the  general  con¬ 
clusions  are  not  affected. 

The  M  factors  covering  representative  ranges  of  h^,  h^,  and  y^,  computed 
from  equations  by  Appleton  and  Beynon,  ^  are  shown  in  Figure  6.  For  the  ratio 
y^/h^,  we  have  chosen  a  range  from  0.4  to  0.6  covering  a  large  range  of  iono¬ 
spheric  conditions.  With  ~  0.4,  the  figure  shows  four  curves  for  h^ 

=  250,  300,  350,  and  400  km.  The  layer  parameters  h^,  h^,  and  y^  used  in  the 
mathematical  formulation  are  interconnected,  y^  and  h^  by  the  selected  ratio, 

and  h  by  the  relation  h  =  h  +  y„.  For  y^/h^  *  0.  6,  the  resulting  M  factor 
m  mo'^m  ■'mo 

5.  Appleton,  E.  V. ,  and  Beynon,  W.  J.G.  (1940)  The  application  of  ionospheric 

data  to  radio  communication  problems.  Part  I, ,  Proc.  Phys.  Soc^  ,  Part  I, 
52:518. 
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HOURLY  CHANGE  IN  MEDIAN  F0F2  (%) 

Figure  5.  Percent  Increase  in  foF2  During  Morning  Hours  for  Low  Sunspot  Activity 


Figure  6.  Dependence  of  M  Factor  on  Skip  Distance  for  '  0*  4  and  0.6 

Respectively 


values  as  a  function  of  range  are  about  2  percent  lower  than  those  for  ~  0*  4. 

Therefore,  only  one  set  of  curves  was  plotted.  The  and  values  computed 

for  each  selected  h_  are  tabulated  separately  for  y_/h  =  0.  4  and  0.  6  in  the  fig- 
m  'mo 

ure.  The  planned  location  of  the  barrier  at  2000-2500 km  for  the  sunrise  period  is 
shown  by  vertical  lines  in  the  figure. 

The  ratios  0.4  and  0.  6  refer  to  thinner  and  thicker  F-layers  respectively. 

With  the  diurnal  variation  of  h^,  higher  at  night  and  lower  at  noon,  the  lower 
curves  represent  night  conditions,  and  the  upper  curves  represent  noon  conditions. 
In  seasonal  variation,  h  is  higher  in  summer  than  in  winter.  Similarly,  with  re¬ 
spect  to  solar  cycle  activity,  the  upper  curves  refer  to  low  sunspot  activity,  and 
the  lower  curves  refer  to  high  sunspot  activity  (at  low  sunspot  activity,  h^  is 
lower  than  al  high  sunspot  activity).  In  descending  order  of  importance,  (he  terms 

dominating  the  M  factor  are  skip  distance,  h  ,  and  y  /h  .  The  slope  of  the  M 

m  m  o 

factor  is  larger  at  small  skip  distances  and  becomes  smaller  as  the  distance  in- 
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creases.  In  Figures  2  through  5  and  the  related  discussions,  we  have  considered 
only  the  variations  of  the  foF2  median  values  at  sunrise  periods  and  have  not  taken 
into  account  the  parameters  h^,  h^,  and  y^,  which  also  change  during  the  sunrise 
period  of  2-3  hr.  During  the  sunrise  period  of  2-3  hr,  the  F-layer  comes  down  by 
about  40  km  (computed  from  lONCAP  model)  at  the  midpoint  of  the  barrier  distance 
(2000-2500  km).  The  M  factor  would  correspondingly  increase  by  10  percent  as 
seen  from  Figure  6. 

For  practical  purposes,  =  0.4  to  0.  6  covers  most  of  the  full  range  of 

ionospheric  changes  expected  during  sunrise  with  four  curves.  One  would  mostly 
be  working  within  the  bounds  of  a  neighboring  pair  of  curves.  For  example,  on  a 
morning  in  a  winter  month  with  high  sunspot  activity,  3SN  *  160,  lONCAP  model 
computations  provide  tlic  following  layer  parameters; 


LMT 

h  (km) 
m 

y^  (km) 
m 

h^  (km) 

y_/h 
m  o 

06 

355 

99 

256 

0.39 

10 

299 

104 

195 

0.53 

To  interpret  these  changes,  one  would  work  in  the  bounds  of  the  middle  pair  of 
curves,  with  the  starting  and  ending  M  factors  of  2.  2  and  2.  9  for  the  midpoint  of 
the  barrier. 

One  way  of  looking  at  the  effects  of  the  changing  ionosphere  on  the  radar  situ¬ 
ation  is  to  start  from  the  relation 


f 


RADAR 

foF2 


M  (Di>i‘ance,  h  ,  h  ,  y  ). 

m  o  m 


(3) 


Keeping  all  layer  parameters  (except  foF2)  and  the  radar  frequency  fixed  and  in¬ 
creasing  (decreasing)  foF2  at  the  midpoint  at  sunrise  (sunset)  causes  the  M  factor 
to  decrease  (increase).  As  Figure  6  shows,  this  causes  the  skip  distance  and 
therefore  the  barrier  to  move  in  (out).  This  fact  combined  with  the  speed  with 
which  foF2  changes  during  sunrise  is  the  main  reason  for  the  need  to  develop  a 
short  term  prediction  scheme. 

To  look  at  the  effect  of  h  ,  let  us  keep  foF2  and  the  radar  frequency,  and 
therefore  the  M  factor,  constant.  Lowering  (raising)  the  layer  height  during  sun¬ 
rise  (sunset)  decreases  (increase  the  skip  distance  and  therefore  the  distance  to 
the  barrier.  Therefore,  the  changes  in  the  F-layer  during  sunrise  (foF2,  h^) 
both  act  in  the  same  direction,  speeding  up  the  movement  of  the  barrier  towards 
the  radar  and  increasing  the  need  to  change  frequency  to  maintain  the  barrier  at 
the  desired  range,  while  providing  overlap  of  part  of  the  radar  footprint  prior  to 
and  after  a  frequency  change. 

From  Figure  6,  the  change  in  M  factor  per  100  km  range  change  is  computed 
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as  a  function  of  skip  distance  and  is  shown  in  Figure  7.  The  figure  shows  that  the 
percentage  change  in  M  factor  as  a  function  of  skip  distance  is  the  same  for  any 
layer  height  for  skip  distances  greater  than  2000  km.  Naturally,  this  does  not 
imply  that  the  M  factor  itself  is  not  changing  as  the  layer  altitude  Increases.  Fig¬ 
ure  6  clearly  shows  the  change  of  M  with  the  change  in  h^. 

The  change  of  M  with  skip  distance  decreases  as  the  skip  distance  increases. 
At  the  center  of  the  barrier,  the  change  of  M  is  2.35  percent  per  100km  move¬ 
ment  of  the  barrier.  Thus,  the  250km  (or  1/2  barrier  width)  allowable  sliding  of 
the  barrier  corresponds  to  a  6  percent  change  in  M  factor.  Therefore,  to  main¬ 
tain  a  50  percent  overlap  of  the  barriers  through  a  frequency  change,  the  radar 
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operating  frequency  has  to  be  changed  whenever  the  change  in  M  factor  is  larger 
than  6  percent.  Figures  5a  and  5b  show  that  the  frequency  has  to  be  changed  5-10 
times  per  hour  to  cover  the  30-60  percent  hourly  change  in  foF2  during  sunrise 
hours.  This  is  equivalent  to  a  change  in  radar  frequency  every  6-12  min.  Addi¬ 
tional  changes  are  required  due  to  the  simultaneous  lowering  of  the  layer.  Such 
frequent  changes  demand  reliable  predictions  of  foF2. 

The  data  base  used  in  the  analysis  is  presented  in  the  following  section. 


.1.  I)  \  l'  \  MASK 

While  defining  the  sunrise  foF2  problem,  we  have  used  foF2  median  values 
from  several  stations  (Figures  2  and  3).  For  investigating  hourly  foF2  prediction 
schemes,  hourly  foF2  data  are  used  from  five  stations:  Narsarssuaq  (NQ),  St. 
Johns  (SJ),  Ottawa  (OT),  Wallops  (WA),  and  Boulder  (BD).  These  data  are  repre¬ 
sentative  of  the  ionosphere  in  the  complete  sector  of  the  East  Coast  OTH-B  cover¬ 
age  area  (see  Figure  1).  The  geographic  and  corrected  geomagnetic  coordinates 
of  all  these  stations  (including  those  in  Figures  2  and  3)  are  listed  in  Table  1. 

Since  foF2  observations  are  routinely  presented  in  UT,  but  the  phenomena  depend 
on  local  time,  the  time  differences  are  shown  in  Table  1.  For  reference,  the 

4 

European  stations  used  by  Davis  and  Rush  for  testing  the  longitudinal  extrapola¬ 
tion  method  are  also  listed  in  the  table.  The  selection  of  stations  for  this  study  is 
guided  by  the  availability  of  hourly  data  and  by  relevance  to  the  ECRS  coverage 
area. 

The  table  shows  that  the  European  stations  are  geomagnetically  further  south 
(equatorwnrd)  than  the  OTIl-I^  coverage  high  latitude  stations.  Since  the  iono¬ 
sphere  is  controlled  by  geogi  aphic  latitude  dependent  phenomena  (solar  production, 
neutral  wind  systems),  as  well  as  by  geomagnetic  latitude  related  phenomena  (K- 
trough,  auroral  oval),  it  is  not  surprising  that  attempts  at  foF2  predictions  for  the 
East  Coast  Radar,  based  on  foE2  tneasurements  in  Europe,  were  not  successful 

A 

(Davis  and  Rush  ).  A  look  at  Table  1  shows  that  the  two  station  method  (Davis  and 
Rush  )  must  extrapolate  foF2  predictions  not  only  beyond  distances  greater  than 
1000  km  but  also  over  time  differences  of  4  hr  for  the  ECRS  application. 

For  the  OTH-relevant  stations,  hourly  foF2  data  are  used  for  four  months: 
March,  June,  September,  and  December,  representing  the  four  seasons.  To  cover 
the  high  and  low  phases  of  solar  cycle  activity,  the  data  for  calendar  years  1969 
(average  smoothed  monthly  sunspot  number  was  100)  and  1975  (average  smoothed 
monthly  sunspot  number  was  15)  arc  used.  Thus,  for  each  station,  the  data  base 
consists  of  8  mo  of  hourly  foF2  observations.  For  prediction,  it  is  essential  to 
have  a  good  data  base.  The  spnrseness  of  the  present  data  base  is  shown  by  tabu- 
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lating  those  periods  when  observations  fall  below  10  in  any  given  month  (Table  2). 

The  table  shows  that  data  suffer  severely  at  Narsarssuaq  and  partially  at  St, 

Johns  and  at  Ottawa.  Most  of  these  data  are  lost  because  parameters  cannot  be 

extracted  from  ionograms  under  spread  F  conditions.  It  is  expected  that  digital 
6  7 

ionosondes  ’  to  be  deployed  for  real  time  OTH  support  in  the  near  future  will 
provide  a  considerably  improved  capability  in  terms  of  spread  F  performance  and 
better  support  to  proposed  prediction  schemes  as  well  as  the  real  time  OTH-B 
operation.  We  do,  however,  feel  that  the  data  base  available  was  sufficient  to  test 
the  various  prediction  schemes. 

The  cumulative  distribution  of  Kp  indices  from  the  respective  4  months  for 
12  hr  around  sunrise  is  shown  in  Figure  8.  The  magnetic  activity  is  not  significantly 


Table  2.  Months  and  Stations  With  Sparse  Data  Base 


1  -  . . .  "  1 

Number  of  Observations 

L.T. 

Station 

Month 

Year 

0123456789  10  11 

1)  Dourbes 

Dec 

69 

9  9 

2)  Narsarssuaq 

Mar 

69 

8  8  9 

Sept 

69 

7 

Dec 

69 

7  5  9 

Mar 

75 

4  2  2  3  6  3  8 

June 

75 

7  3 

Sept 

73 

10  13  3  7 

Dec 

75 

22467864 

3)  St  Johns 

Mar 

75 

7  8  5  4 

Dec 

75 

4  6  7  7  5 

4)  Ottawa 

Mar 

75 

9  9  7  5  G 

Dec 

75 
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Kp  INDEX 


Figure  8.  Cumulative  Distribution  of  Kp  Index  for 
the  Period  Selected  for  Study 


different  during  the  periods  of  high  and  low  solar  cycle  activity  selected  in  this 

O 

study.  Mendillo  et  al  showed  that  the  effect  of  magnetic  activity  on  foF2 
is  not  clearly  understood  and  possibly  depends  on  the  start  time  (LT)  of  th,*  mag¬ 
netic  disturbance.  They  inspected  28  storms  and  found,  that  after  a  short  enhance- 


8.  Mendillo,  M. ,  Lynch,  F.X.,  and  Klobuchar,  J,  A.  (1980)  Geomagnetic  activity 
control  of  ionospheric  variability,  in  Prediction  of  Terrestrial  Effects  of 
Solar  Activity,  Vol.  4  of  Solar-Terrestrial  Predictions  Proceedings,  R.F, 
Donnelly,  Ed.,  Space  Environment  Laboratory,  National  Oceanic  and  Atmos¬ 
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ment  at  the  start  of  the  storm,  foF2  is  subsequently  depressed  over  a  period  of 
several  days.  The  present  study  confirms  these  findings,  showing  that  on  highly 
disturbed  days  [5  per  month  identified  by  World  Data  Center  (WDC),  Boulder,  Co¬ 
lorado],  the  respective  foF2  values  are  usually  below  the  observed  monthly  median 
foF2s.  However,  a  quantitative  relation  between  Kp  and  the  respective  hourly 
foF2  has  not  been  found. 

Since  we  plan  to  use  predictions  from  ionospheric  models  in  the  analysis,  the 
sunspot  number  predicted  by  WDC  Boulder  (1  month  in  advance),  and  the  observed 
sunspot  numbers  for  the  8  mo  covered  in  the  analysis,  are  listed  in  Table  3.  Be¬ 
tween  high  and  low  solar  activity,  the  sunspot  numbers  changed  by  a  factor  of  10. 


Table  3.  Sunspot  Activity  for  the  Selected  Period 


Year 

Month 

Prediction 

(1  Month  Before) 

Ob¬ 

served 

1969 

March 

99.  8 

135,8 

June 

106.  3 

106.0 

September 

101.  9 

91.  3 

December 

96.4 

97.9 

1975 

March 

23.0 

11.5 

June 

18.  3 

11.4 

September 

15.  7 

13.9 

December 

14.4 

7.  8 

I.  VWI.^SIS  \M)  DISCI  SSION 

WTiile  considering  the  constraints  in  maintaining  the  barrier  range,  we  have 
seen  that  the  major  concern  is  in  the  rate  of  change  in  foF2,  rather  than  in  the  ab¬ 
solute  value  of  foF2.  The  variation  in  foF2  for  two  stations,  St.  Johns  and  Wal¬ 
lops,  representative  for  the  northern  and  central  ionosphere  of  the  ECRS,  is 
shown  in  Figures  9  to  12.  Similar  variations  are  observed  at  all  the  stations. 

Each  of  these  figures  is  divided  into  two  sections  along  the  vertical.  The  left- 
hand  section  presents  results  from  analysis  for  St.  Johns  and  the  right-hand  sec¬ 
tion  shows  results  for  Wallops. 

The  diurnal  variation  of  foF2  during  the  morning  transition  for  the  equinoxes 
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WflLLOPS 


Figure  9a.  Diurnal  Variation  of  foF2  in  the  Morning  for  Eqinox 


Figure  9b.  Diurnal  Variation  of  foF2  in  the  Morning  for  Solstice 


variability 


ST  JOHNS  WALLOPS 


L  S  T 

Figure  10b.  Percent  Variation  of  foF2  in  the  Morning  Hours  for  Solstice 
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ST  JOHNS 


Figure  11a.  Hourly  Change  (Slope)  of  foF2  in  the  Morning  Hours  for  Equinox 
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Figure  lib.  Hourly  Change  (Slope)  of  foF2  in  the  Morning  Hours  for  Solstice 
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Figure  12a.  Percent  Hourly  Variation  of  foF2  in  the  Morning  Hours  for  Equinox 


of  high  (1969)  and  low  (1975)  sunspot  activity  is  shown  in  Figure  9a  and  for  the  sol¬ 
stices  in  Figure  9b,  The  curves  present  the  median  and  to  values  of  foF2.  A 
feature  common  to  all  the  curves  is  the  rapid  increase  in  foF2  during  the  sunrise 
period.  The  effect  is  largest  for  the  winter  month  of  December  and  smallest  for 
the  summer  month  of  June.  The  absolute  values  of  foF2  for  solar  maximum  (1969) 
are  twice  as  large  of  those  for  solar  minimum  (1975).  During  peri'^ds  of  low  solar 
activity  (1975),  the  foF2  values  at  the  geomagnetically  high  latitude  station  St. 
lohns  are  routinely  2  MHz  or  lower,  before  sunrise.  These  consistently  low  foF2s 
would  barely  allow  the  establishment  of  a  barrier  at  the  nominal  range  of  2000  km 
(1100  nm)  during  the  late  night  period,  at  low  solar  activity.  For  all  the  repre  - 
sentaiivo  months,  there  is  a  rapid  increase  (large  slope)  in  foF2  for  the  2  hours 
at  sunrise  followed  by  .s  relatively  slower  increase  (smaller  slope)  for  2  more 
hours. 

The  relation 


P  = 


100  X 


(foF2+q  -  foF2.g) 
2  foF2  Med 


(4) 


is  used  to  compute  the  percentage  relative  variability  in  foF2  for  each  hour.  Re¬ 
sults  are  shown  in  Figures  10a  and  10b.  The  curves  show  that  the  percent  varia¬ 
tion  does  not  display  any  systematic  diurnal,  seasonal,  or  solar  activity  control. 
The  relative  day-to-day  <  ariability  of  the  hourly  data  is  (  f  the  order  cf  15-20  per¬ 
cent  around  the  median  Va  ues  of  foF2.  This  variability  is  2  or  3  times  larger  than 
that  allowable  (~6  percent)  to  prevent  the  barrier  from  moving  by  more  than  250  km 
(1 '2  of  its  width).  Figures  10a  and  10b  show  that,  using  a  fixed  (for  example,  pre¬ 
dicted  from  lONCAP  for  a  given  hour)  operational  frequency,  the  barrier  may 
start  anywhere  in  at 750 km  distance  from  the  planned  2000km  nominal  distance. 
This  15  percent  variability  around  th*"  hourly  median  foF2  is  superposed  on  the  20- 
50  percent  hourly  change  in  foF2  during  the  sunrise  transition  hours.  It  is,  there¬ 
fore,  essential  to  reduce  this  uncertainty  through  a  short  term  prediction  scheme. 

For  an  effective  real-time  short  term  prediction,  it  is  necessary  to  estimate 
the  expected  hourly  changes  in  foF2.  The  hourly  slope  for  the  two  stations  is  pre¬ 
sented  in  the  same  format  in  Figures  11a  and  11b.  In  addition  to  median  and  d  o 
levels  of  the  hourly  slope  of  foF2,  the  lowest  slope  observed  for  the  month  for  the 
particular  hour  is  shown  in  the  figure.  The  hourly  change  in  foF2  is  largest  in 
winter  months  and  smallest  in  summer  months.  The  change  in  slope  is  a  factor  of 
2  larger  during  the  high  phase  of  solar  activity  than  during  the  lo  v  phase  of  solar 
activity.  The  lowest  slope  is  observed  to  occur  on  the  days  of  magnetic  storms. 
The  highest  slope  (not  shown  in  the  figures)  occurs  on  magnetically  quiet  days 
and/or  on  days  corresponding  to  the  recovery  from  magnetic  storms.  Even  the 
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curve  of  lowest  slope  shows  positive  magnitude  for  a  few  hours  at  sunrise  indicat¬ 
ing  that  the  increase  in  foF2  at  sunrise  is  stronger  than  the  decrease  in  foF2  asso¬ 
ciated  with  magnetic  storms. 

The  relative  (percent)  hourly  slope  is  computed  using  the  relation 


Slope  *  200  X  ■ 


(foF2i  +  1  -  foF2i) 


(foF2£  +  1  +  foF2£) 


(5) 


Plotted  in  Figures  12a  and  12b  are  the  values  of  the  median  observed  slope  and  the 
±0  level  of  the  slope  distribution.  The  figures  show  that  the  relative  slope  is 
larger  in  winter  months  (December)  by  almost  a  factor  of  2,  as  compared  with 
summer  and  equinox  data.  Hut  there  is  no  significant  difference  in  relative  slope 
for  the  respective  seasons  between  high  and  low  solar  activity.  The  relative  slope 
shows  a  strong  diurnal  variation.  The  relative  slope  reaches  its  peak  2  hr  after 
sunrise  and  decreases  slowly  for  the  following  2  hr.  Thus,  the  foF2  changes  asso¬ 
ciated  with  sunrise  are  significant  for  a  duration  of  3-4  hr.  The  typical  peak  me¬ 
dian  slope  of  30  percent  suggests  that  the  radar  frequency  has  to  be  adjusted  on  the 
average  5  times/hour  to  accommodate  the  hourly  change  in  foF2  during  the  sunrise 
period. 

All  the  above  results  can  be  briefly  summarized  as:  (1)  The  foF2  noon  values 
as  well  as  the  hourly  slope  of  foF2  during  sunrise  hours  are  largest  in  winter  and 
smallest  in  summer:  (2)  The  relative  magnitude  (percent  increase)  of  the  hourly 
change  in  foF2  is  twice  as  large  in  winter  as  during  equinox  and  summer;  (3)  The 
relative  magnitude  of  the  hourly  change  in  foF2  for  the  particular  seasons  is  rea¬ 
sonably  independent  of  solar  cycle  activity;  (4)  The  range  of  foF2  variability  at  any 
given  hour  (  >  15-20  percent)  is  larger  than  the  range  acceptable  (<  6  percent)  for 
the  satisfactory  positioning  of  the  barrier;  (5)  This  variation  is  superposed  on  the 
rapid  systematic  increase  in  foF2  (of  >  30  percent  per  hour)  observed  for  about 
4  hr  during  the  sunrise  period,  for  6  mo  covering  the  seasons  of  late  fall-winter 
to  early  spring. 

One  of  the  approaches  to  improving  the  prediction  of  foF2  at  one  station  is  to 
use  and  '’rtrapolate  recent  observations  from  another  station  to  the  east  of  the  lo¬ 
cation  of  interest.  Th*’  method  can  be  called  the  "two  station  method.  "  For  the 
OTH-B  coverage,  European  stations  such  as  Dourbes  (Belgium),  Slough  (England), 
or  Lindau  (Germany)  could  potentially  satisfy  the  requirements  for  such  predic¬ 
tions  at  St.  Johns  in  the  OTH-B  coverage  area,  provided  the  two  station  method 
works.  When  one  compares  the  sunrise  foF2  slopes  over  a  few  hours  from  two 
stations  like  St.  Johns  and  other  European  stations,  one  otitains  a  significant  cor- 
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relation  (  ^0.4),  because  of  the  large  systematic  increase  in  foF2  due  to  the  sun¬ 
rise  effect  at  both  stations. 

3 

Rush  and  Miller  analyzed  foF2  data  from  32  ionospheric  stations  to  deter¬ 
mine  the  correlation  distances  for  foF2  for  the  same  local  time  from  1958  data. 
They  found  that  the  correlation  distance  along  the  E-W  direction  is  larger  than  the 
correlation  distance  in  the  N-S  direction.  They  found  a  diurnal  dependence  of  the 
correlation  distance,  with  largest  correlations  at  06  LMT  and  smallest  at  00  LMT. 
The  correlation  distance  at  06  LMT  was  2200  km  in  the  E-W  direction  and  1500km 
in  the  N-S  direction.  The  separation  between  European  stations  Llndau,  Slough, 
or  Dourbes  and  the  midpoint  sounder  at  St.  Johns  is  greater  than  4000  km.  This 
distance  is  much  larger  than  the  E-W  correlation  distance  of  2200  km  beyond  which 

4 

extrapolation  becomes  unreliable.  Davis  and  Rush  analyzed  the  data  for  the 
European  stations  of  OTH-B  interest  and  confirmed  the  previous  findings  that  on 
the  basis  of  distance  alone,  the  foF2  observations  cannot  be  extrapolated  from 
European  stations  to  the  OTH-D  area  for  obtaining  reliable  predictions. 

For  this  study,  the  period  of  sunrise  is  determined  to  last  for  5  hr  (six  hourly 
observations)  consisting  of  1-2  hr  of  pre-sunrise  characterized  by  a  flat  shape  of 
the  foF2  curve  followed  by  3-4  hr  continuous  increase  in  foF2.  The  foF2  predic¬ 
tions  were  tested  for  all  data  available  for  these  six  hourly  observations. 

In  this  study,  we  have  investigated  the  use  of  real  time  and  recent  history  of 
foF2  data  from  a  set  of  five  ionospheric  stations  (Narssarssuaq,  St.  Johns,  Ottawa, 
Wallops  Island,  and  Boulder)  for  the  short  term  prediction  (1  hr  in  advance)  of  the 
rapid  foF2  changes  during  sunrise.  We  considered  this  approach  because.  If 
found  useful,  it  would  allow  the  prediction  of  the  foF2  sunrise  transition  using  the 
chain  of  ionospheric  sounders  deployed  for  real  time  ionospheric  specification  of 
the  ECRS  midpoint  ionosphere  (Goose  Bay,  St.  Johns,  Bermuda,  Wallops  Island). 

To  judge  the  performance  of  the  various  prediction  schemes,  we  established 
as  a  baseline  the  prediction  errors  that  resulted  from  the  comparison  of  actual 
hourly  measurements  with  the  predicted  foF2  determined  from  the  lONCAP  model 
using  the  predicted  (1  month  in  advance,  not  actual)  monthly  smoothed  sunspot 
number  for  the  respective  month. 

The  proposed  prediction  scheme  uses  recent  history  (for  past  n  days)  and  real 

2 

time  foF2  observations  (day  n  +  1).  Various  combinations  of  lONCAP  model  pre¬ 
dictions  and  foF2  observations  are  also  tested  to  determine  the  best  possible  pre¬ 
diction  to  cover  the  case  where  no  real  time  data  are  available  due  to  interruptions 
in  the  ionosonde  operations. 

The  proposed  scheme  is  based  on  the  use  of  the  foF2  empirical  data  base  and 
is  shown  schematically  in  Figure  13,  The  prediction  uses  averaged  foF2  data  for 
four  preceding  days.  As  a  first-order  correction  the  predicted  foF2  is  equal  to 
the  foF2  now  (observed)  plus  the  average  foF2  slope  for  the  next  hour.  As  a 
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1st  ORDER  foFa  NOW + 1  HR  SLP. 

2nd  ORDER  foF2  NOW+IHR  SLP.  +  SLP.  CORR. 


Figure  13.  Scheme  for  Prediction  of  foF2  at  Hourly  Intervals  From  Observed 
"Now  foF2"  and  Hourly  Slope  of  foF2 


second-order  correction  a  slope  correction  term  is  made  proportional  to  the  dif¬ 
ference  between  foF2  now  and  the  averaged  foF2  (for  four  previous  days). 

We  will  first  discuss  the  mathematical  procedures  used  for  the  various  ap¬ 
proaches  tested  and  subsequently  we  will  discuss  the  performance  of  each  scheme. 

The  prediction  for  the  day  n  +  1  (prediction  day)  based  on  the  averaging  of 
samples  from  n  preceding  days  for  the  corresponding  hour  i  is  given  by 

n 

_ _  E  foF2n.  1 

.  1.  I'  ■  i  •  ■‘—5 - 

with  p(loF2)  the  predicted  foF2  value.  Kq.  (6)  provides  averaged  foF2  for  the  past 
n  days  as  the  zeroth  order  prediction. 

If  the  past  days'  data  is  not  available,  the  prediction  has  to  be  made  using  the 
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lONCAP  model,  so  that  the  right-hand  side  of  Eq.  (6)  is  simply  the  foF2  predicted 
for  that  hour  by  the  model. 

Since  real  time  hourly  (or  half  hourly)  foF2  data  will  be  available  to  the  En¬ 
vironmental  Assessment  (EA)  operator  at  the  radar  site,  we  tested  the  concept  of 
using  the  actual  (i  hour)  measurement  of  foF2  in  combination  with  an  average  slope 
derived  from  the  prior  days'  average  of  foF2  measurements  at  the  i  and  i  +  1  hour 
to  predict  foF2  for  the  i  +  1  hour  for  day  n  +  1  (prediction  day).  This  scheme  is 
expressed  formally  as 


p(fol'2 


1.  i 


+  1, 


.  + 
1 


n 

E  (foF2 

1 


n*  1  + 


1  -  j) 


(7) 


with  n  -  1  as  tlie  day  for  which  the  prediction  will  be  made,  and  n  as  the  number  of 
previous  days  from  which  the  average  foF2  has  been  determined  for  the  ith  hour. 

In  case  of  lack  of  the  past  data,  the  second  term  on  the  right-hand  side  of  Eq.  (7) 
can  be  replaced  by  the  lONCAP  slope. 

Finally  we  added  a  second  order  slope  correction,  which  depends  on  the  mag¬ 
nitude  of  the  current  foF2  in  relation  to  the  past  n  day  average  foF2.  The  predic¬ 
tion  is  given  by; 


n 

I  (foF2„,  i  ^  1  -  foF2n,  j)  ^ 

1,  i-  l’  1,  i^  " 

^  i  -  foF2n  +  l.J,  -j  (gj 

^0^2^,  i 

This  means  that  when  the  observed  foF2  is  below  the  average  foF2  of  the  past 
n  days,  the  slope  correction  is  larger  than  the  average  hourly  slope.  If  foF2  is 
above  the  average  foF2,  the  correction  is  smaller  than  the  average  hourly  slope. 
This  correction  assumes  that  the  tendency  is  to  return  towards  the  median/average 
curve  (value). 

Note  that  the  averaged  term  in  the  square  bracket  is  same  as  that  in  Eq.  (6). 

In  Eqs.  (6)-(8),  one  can  use  foF2  observations  and/or  foF2  predictions  from 

the  lONCAP  model.  In  general,  the  analysis  shows  that  the  use  of  recent  and  real 

time  observations  in  Eqs.  (6)-(8)  provides  predictions  better  than  those  obtained 

2 

by  using  foF2s  from  the  lONCAP  model. 
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The  method  of  slopes  of  foF2  was  used  for  predictions  at  intervals  of  1,  2, 
and  3  hr  respectively.  As  was  expected,  the  analysis  using  the  1-hr  interval  re¬ 
sulted  in  the  smallest  errors.  Only  the  results  for  a  1-hr  interval  are  considered 
in  the  following  summary.  A  slope  correction  similar  to  Eq.  (8)  to  account  for  the 
level  of  magnetic  activity  [foF2  slope  “  (Kp^  ^  j-Kp)]  was  tested.  (The  WDC  Kp 
index  is  available  at  3-hr  intervals  only. )  The  correction  was  found  to  be  insigni¬ 
ficant.  Therefore,  no  such  correction  was  used  in  the  prediction  scheme. 

The  average  of  prior  n  days  as  a  prediction  was  used  with  values  of  n*  1,  2, 

3,  4,  5,  10,  15,  20,  25.  It  was  determined  that  for  n  =  4  the  prediction  of  foF2 
yields  the  least  error.  Therefore  in  the  following  discussion  only  n  =  4  is  used. 

Eqs.  (6),  (7),  and  (8)  provide  a  step  by  step  approach  for  improving  foF2  pre¬ 
dictions.  The  improvement  of  one  prediction  scheme  over  the  other  is  checked  in 
terms  of  reduction  in  relative  error  [see  Eq.  (9)],  reduction  in  AfoF2  [see 
Eq.  (10)],  and  increase  in  correlation  coefficient  showing  higher  percent  of  cor¬ 
related  population  [see  Eq.  (12)].  First  we  will  consider  the  relative  error  in 
prediction. 

The  improvement  by  a  given  prediction  scheme  is  determined  by  comparing 
errors  in  prediction  at  the  median  (50  percent  level),  ±a  level  and  the  total  range 
of  the  data  population.  For  each  observation  the  percent  error  is  determined  by 


foF2obs  “  ^°^2p|.e(jicted 

4foF2  (%)  =  - - - 

foF2obs 


(9) 


The  errors  are  computed  from  the  hourly  observations,  covering  a  6-hr  period 
(2  hr  before  and  4  hr  after  sunrise),  for  each  month.  From  this  monthly  set  of 
error  data,  values  are  determined  for  median,  125  percent,  ±o  ,  ±1.5o,  12  o  levels 
and  total  coverage  (extreme  values)  of  the  monthly  data.  In  the  tables,  3  levels 
are  listed  for  discussion:  median, io  ,  and  total  coverage  range. 

As  an  example,  the  results  for  Boulder  for  June  1969  are  shown  graphically 
in  Figure  14.  The  figure  includes  prediction  errors  to  cover  ±25  percent,  to  , 

41.  5o,  and  i2o  levels  of  the  error  populations.  In  the  figure  for  any  given  level  of 
the  populations,  the  errors  are  largest  for  lONCAP  predictions  and  smallest  for 
the  proposed  method  using  second-order  corrections  [see  Eq.  (8)].  The  median 
error  from  iONCAP  is  biased,  indicating  the  inability  of  lONCAP  to  predict  the 
median  foF2  for  the  sunrise  period,  even  if  the  sunspot  number  is  predicted  cor¬ 
rectly  (June  1969,  see  Table  3).  The  errors  are  progressively  reduced  as  one 
uses  the  past  4  days'  (hourly)  averages  of  foF2,  the  foF2  now  +4  days'  average 
slope,  and  adding  the  modified  slope  [liq.  (8)].  An  important  aspect  of  the  use  of 
the  recent  past  and  the  current  measurement  is  that  the  error  distribution  is  very 
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Figure  14.  Comparison  in  the  Prediction  Error  (Percent)  of  the  foF2  for  Various 
Prediction  Schemes 


nearly  centered  on  zero,  that  is,  in  the  median  over  a  month,  the  prediction  is 
close  to  perfect.  Therefore,  these  approaches  make  the  prediction  task  indepen¬ 
dent  of  specifying  or  deriving  a  sunspot  number,  to  drive  a  prediction  model. 

The  comparison  of  the  percent  error  of  the  various  methods  with  the  scatter 
of  the  data  base  used  around  the  data  median  (see  last  column  of  Figure  14)  is 
identical  to  a  comparison  with  the  hypothetical  case,  that  lONCAP  could  predict 
the  monthly  median  (for  the  sunrise  hours)  exactly.  The  "improvement"  of  the 
proposed  scheme  over  the  scatter  distribution  of  the  data  base  is  also  a  good  esti¬ 
mate  of  the  improvement  over  the  hypothetical  case  of  correct  median  prediction 
l)y  lONCAP.  The  comparison  shows  that  the  spread  of  the  observational  data 
around  the  median  is  larger  than  the  spread  of  the  errors  from  the  proposed  pre¬ 
diction  scheme  (Figure  14).  The  figure  emphasizes  the  fact  that  using  the  observed 
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foF2  in  Eqs,  (6)-(8)  yields  better  predictions  than  those  based  on  the  lONCAP 
(ITS-78)  model.  ^ 

The  results  shown  in  Figure  14  generally  repeated  in  the  evaluation  of  the 
five  station/8  months  study.  We,  therefore,  select  Eq.  (8)  as  the  proposed  scheme 
for  improved  foF2  prediction  during  sunrise.  In  the  subsequent  discussion  of  the 
prediction  improvement,  we  will  limit  ourselves  to  the  comparison  of  the  errors 
resulting  from  the  proposed  method  with  those  from  lONCAP  and  with  the  scatter 
of  the  data  base  around  the  median.  This  will  eliminate  the  listing  of  all  results 
showing  performance  worse  than  that  by  Eq.  (8). 

The  results  from  the  analysis  are  summarized  in  Table  4,  The  first  column 
lists  the  level  of  the  population  at  which  the  errors  are  computed.  The  second  col¬ 
umn  lists  the  prediction  scheme.  The  third  "method"  is  essentially  the  spread  in 
the  observed  data  base.  The  third  column  identifies  the  ionospheric  stations  pro¬ 
viding  the  data.  The  next  eight  columns  present  the  relative  (percent)  error  in 
prediction  for  four  seasons  for  high  and  low  sunspot  activity  periods.  The  last 
column  presents  the  average  error  value  for  a  given  station.  Looking  across  the 
columns,  no  seasonal  dependence  is  evident,  suggesting  the  extent  of  data  spread 
(relative)  around  the  median  value  is  independent  of  seasonal  and  solar  (sunspot) 
variation.  The  comparison  of  average  errors  from  the  last  column  for  the  predic¬ 
tion  schemes  shows  that  errors  are  significantly  smaller  for  the  proposed  method 
at  all  levels  of  the  distribution. 

A  reduction  in  error  by  the  proposed  method  from  those  from  lONCAP/ob- 
served  data  set  is  computed  by 

4 error  *  | ERROR  from  lONCAP j  -  | ERROR  from  proposed  method]  (10) 

at  respective  levels  (from  Table  4).  The  results  are  summarized  in  Table  5.  In 
the  table,  a  negative  sign  preceding  the  value  shows  that  the  suggested  method 
caused  a  deterioration  instead  of  Improvement  over  the  method  under  consideration 
(blocks  with  odd  sequence  number  refer  to  lONCAP,  even  number  to  observational 
data  set).  A  look  at  the  last  column  shows  that  the  average  prediction  error  is  re¬ 
duced  by  8,  8,  and  13  percent  over  lONCAP  at  the  levels  of  the  median,  ±o  ,  and 
total  population.  In  the  observed  data  set,  the  first  term  on  the  right-hand  side  of 
Eq.  (10)  is  the  spread  in  data  at  the  respective  levels  of  population.  Over  the  ob¬ 
served  data  set,  the  prediction  error  is  reduced  by  4  and  7  percent  in  comparison 
to  the  natural  spread  at  ±la  and  total  population,  respectively.  The  result  shows 
that  predictions  by  the  proposed  scheme  are  better  than  the  scatter  of  the  observed 
data  (that  is,  predictions  by  a  perfect  lONCAP  model). 

The  percentage  reduction  in  error  by  the  proposed  method  with  respect  to 
lONCAP  and  the  observed  data  set  can  be  computed  by 
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Table  4.  Error  (Percent)  in  foF2  for  Various  Prediction  Methods 
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Table  4  (coiit).  Krror  (Percent)  in  foF2  for  Various  Prediction  Methods 
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Table  4  (cont).  Krror  (Percent)  in  foF2  for  Various  Prediction  Methods 
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Table  5.  Error  Reduction  (Percent)  by  Proposed  Method  Over  Other  Methods 
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(4ERRORi  -  4ERR0R proposed  METHOD^ 

factor  (percent)  =  100  x - - 

ERROR  [ 


(11) 


where  i  -  refers  to  values  from  lONCAP  and  from  observational  data  set  respec¬ 
tively.  The  term  dERROR  is  the  spread  in  the  observational  data  base. 

The  improvement  factors  are  computed  by  use  of  respective  values  from 
Table  4.  The  results  from  the  last  columns  of  Tables  4  and  5  and  evaluation  by 
Eq.  (11)  arc  summarized  in  Table  6.  The  table  is  divided  into  three  major  sec¬ 
tions.  These  list  error,  reduction  in  error  and  improvement  respectively.  Each 
column  is  subdivided  in  three  columns  to  present  the  median,  ±1  o  and  total  pop¬ 
ulation  coverage.  Each  block  (along  the  row)  refers  to  a  prediction  method.  The 
line  at  the  bottom  of  each  block  presents  the  average  value  from  all  the  stations. 
Referring  to  lines  in  the  last  three  columns,  the  proposed  method  shows  86  per¬ 
cent  improvement  (factor)  in  prediction  over  that  from  lONCAP.  Such  a  large  im¬ 
provement  is  due  to  a  reduction  in  bias  in  the  lONCAP  prediction  seen  in  Figure  14. 
At  the  il  0  level,  the  proposed  method  provides  an  improvement  (factor)  of  47  per¬ 
cent  ov'er  lONCAP  and  of  29  percent  over  the  observed  data  set.  The  overall  (to¬ 
tal  population)  improvement  (factor)  is  26  percent  over  lONCAP  and  of  16  percent 
over  the  observed  data  set. 

One  may  choose  the  reduction  in  AfoF2  instead  of  that  relative  to  foF2  (as  pre¬ 
sented  in  Tables  4  to  6).  The  results  are  presented  in  Tables  7  to  9  respectively, 
using  the  same  format  as  for  Tables  4  to  6, 

A  look  at  the  first  block  of  Table  9,  presenting  the  summary  of  the  results, 
shows  that,  on  the  average,  lONCAP  predictions  are  off  by  iO.  4,  ±  0.  7,  and  ±1.7 
MHz  at  median,  ilo  and  total  population  levels.  The  improvement  by  the  sug¬ 
gested  method  over  lONCAP  reduces  the  error  at  these  levels  by  u.  4  MHz.  The 
last  block  shows  that  the  prediction  by  the  new  method  reduces  the  error  compared 
to  the  spread  in  the  natural  data  base  by  0.  2  Mllz  at  the  tlo  level  and  by  0.3  MHz 
in  the  overall  data  base. 

The  second  block  shows  that  the  new  method  reduces  the  prediction  error  to 
±0.  4  MHz  at  the  ±o  level  population  satisfying  the  prediction  requirement  of  ±0.  5 
MHz  suggested  for  the  East  and  West  Coast  radar  operations. 

The  tables  show  that  the  method  proposed  here  will  provide  better  predictions 
than  lONCAP,  and  also  a  lesser  scatter  in  prediction  errors  than  those  existing  in 
the  natural  data  base  (that  is,  improved  prediction  over  that  from  a  perfect 
lONCAPi. 

The  analysis  of  the  various  prediction  schemes  was  extended  to  restrict  pre¬ 
diction  to  quiel  days  only.  The  data  of  the  five  most  disturbed  days  in  each  month 
were  excluded  from  the  analysis,  and  no  predictions  were  computed  for  the  dis- 


46 


Table  6.  Summary  of  Error,  Reduction  in  Error  and  Improvement  by  Proposed  Method 
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Table  8.  Reduction  in  dfoF2  (MIlz)  by  Proposed  Method  Over  Other  Methods 
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Table  9.  Summary  of  /lfoF2  (MIlz)  Error  and  Reduction  in  Error  (dfoF2)  by 
Proposed  Method 


d  foF2  MHz  at  Level 

dfoF2  MHz  (Reduction)  at  Level 

Method 

St.  No 

El 

Extremes 

Median 

±0 

Extremes 

lONCAP 

1 

0.4 

0.6 

1.5 

0.3 

0.2 

0.5 

0.3 

0.8 

1.9 

0.2 

0.3 

0.5 

0.4 

0.8 

1.9 

0.4 

0.4 

0.4 

0.6 

1.5 

0.2 

0.2 

0.3 

5 

0.7 

0.9 

1.9 

0.6 

0.5 

0.5 

Av 

0.4 

0.7 

1.7 

0.4 

0.3 

0.4 

Proposed 

1 

0.1 

0.4 

1.0 

2 

0 

0.4 

1.4 

3 

0 

0.4 

1.5 

4 

0.1 

0.4 

1.2 

5 

0 

0.4 

1.4 

Av 

0.1 

0.4 

1.3 

Observed 

1 

0 

0.5 

1.4 

0.1 

0.4 

Data  Set 

2 

0.1 

0.6 

1.8 

0.  1 

0.3 

3 

0 

0.7 

1.8 

0.3 

0.3 

4 

0 

0.5 

1.4 

0.1 

0.2 

5 
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0.6 

1.8 

0.2 

0.3 

Av 

0 

0.6 

1.5 

0.2 

0.3 

1.  Narsarssuaq  2, 

St.  Johns  3. 

Ottawa  4. 

Wallops 

5.  Boulder 

turbed  days'  foF2s.  The  results  are  slightly  better  but  not  significantly  different 
from  those  presented  in  Tables  4  to  9, 

In  the  above  discussion,  v/e  considered  relative  (percent)  and  magnitude  im¬ 
provements  in  foF2,  obtained  by  the  proposed  method  over  the  lONCAP  method 
and  over  the  observational  data  set. 

4 

The  degree  of  improvement  equation  used  by  Davis  and  Rush  deals  with  the 
fractions  of  correlated  and  uncorrelated  populations  of  the  data  set.  The  equation 

9 

(Gautier  and  Zacharisen  )  is 


9.  Gautier,  T.N.,  and  Zacharisen.  D.H.  (1965)  Use  of  space  and  time  correla¬ 
tions  in  short  term  ionospheric  predictions.  Conference  Record  1st  Annual 
IEEE  Communications  Convention,  pp.  671-676. 
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I  =  100  [l  -  (1  -  R^) 


(12) 


where  I  is  the  degree  of  improvement  in  percent  and  R  is  the  correlation  coeffi¬ 
cient. 

First  let  us  delete  the  multiplying  factor  100  so  that  R  and  I  become  compar¬ 
able  quantities  with  ranges  from  -1  to  +  1.  Further  replacing  I  by  1  -  I',  we  obtain 
2  2 

R  +  r  =1  where  the  partition  of  the  correlated  and  uncorrelated  fractions  of  the 
data  set  is  clearly  evident. 

The  distributions  of  these  correlated  and  uncorrelated  populations  for  the  two 

prediction  schemes:  (1)  lONCAP,  (2)  Proposed  Method,  and  for  the  use  of  the 

scatter  of  the  observed  data  around  the  monthly  median  foF2  (perfect  lONCAP)  are 

computed.  The  least  squared  deviation  with  first  order  (linear)  fit  is  computed 

among  foF2s  predicted  by  the  proposed  method  vs  the  observed  hourly  foF2  data 

set.  The  correlation  coefficients  are  summarized  in  Table  10,  which  has  the 

same  format  as  that  for  Table  9.  The  third  column  presents  the  values  for  the 

correlation  coefficients.  The  range  of  percent  population  correlated  is  given  by 
2 

R  to  R  as  listed  in  the  next  column.  The  table  shows  that  the  correlation  coeffi¬ 
cients  for  the  proposed  method  are  highest  compared  to  those  from  the  lONCAP 
and  the  observed  median  foF2s.  The  degree  of  improvement  computed  by  the  use 
of  Eq.  (7)  is  in  the  following  column.  These  again  confirm  that  the  proposed 
method  offers  an  Improvement.  The  improvement  (dl)  achieved  by  proposed 
method  over  lONCAP  and  the  observed  data  distribution  (perfect  lONCAP)  is  pre¬ 
sented  in  column  6.  The  results  show  that  additional  14  and  10  percent  populations 
are  shifted  from  the  non-correlated  group  to  the  correlated  group,  for  lONCAP 
and  observed  median  foF2  prediction  methods.  The  last  column  shows  relative 
improvement  [similar  to  Eq.  (11)].  Thus  the  proposed  method  offers  a  relative 
improvement  of  29  and  20  percent  over  the  respective  methods.  The  transmitter 
site  is  1100-1340  km  from  the  mid-point  region  of  OTH-B  interest.  Though  the  cor¬ 
relation  distance  is  just  at)out  right,  tlie  two  station  method  is  unable  to  predict 
hourly  foF2  at  St.  Johns  for  the  sunrise  transition  period  from  the  larger  and  lead¬ 
ing  changes  in  foF2  observed  at  the  European  stations. 


CONCI.I  SIONS 

The  analysis  of  data  from  five  stations  representing  the  OTH-B  coverage  area 
has  been  conducted  for  four  seasons  of  two  selected  years  to  cover  high  and  low 
phases  of  solar  cycle  activity. 

Between  two  prediction  schemes:  (1)  use  of  data  from  a  single  station,  and 
(2)  extrapolation  of  data  from  one  station  to  another,  the  first  is  better.  For  a 
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Table  10.  Comparisons  of  Correlations  and  Improvements  Among  Prediction  Schemes  Used 
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Narsarssuaq  2.  St.  Johns  3.  Ottawa  4.  Wallops  5.  Boulder 


given  station,  out  of  various  alternatives,  such  as  the  use  of  observed  median 
foF2,  lONCAP  predictions  using  the  predicted  sunspot  number,  and  averaging  the 
data  for  "n"  number  of  days  prior  to  the  prediction  day,  the  four-day  averages  is 
the  best  approach.  The  prediction  is  further  improved  by  application  of  first-  and 
second-order  slope  correction  shown  by  Eqs.  (7)  and  (8).  The  predictions  are 
very  slightly  improved  if  magnetically  disturbed  days  are  neglected  in  the  predic¬ 
tion  scheme.  An  important  feature  of  the  proposed  scheme  is  that  the  prediction 
error  distribution  has  a  negligible  median  error.  The  scheme  makes  short  term 
predictions  independent  of  the  still  imprecise  prediction  of  the  relevant  sunspot 
numbers. 

In  summary,  the  steps  in  the  prediction  of  foF2  during  the  sunrise  period  are 
as  follows:  (1)  For  each  LT  hour  of  the  sunrise  period,  store  the  hourly  foF2s 
from  observations  of  up  to  4  days  prior  to  prediction.  Exclude  the  days  that  are 
found  to  be  magnetically  severely  disturbed.  (2)  Determine  the  hourly  slope  of 
foF2  for  each  day  and  average  up  to  four  slopes.  (3)  To  the  current  foF2  obser¬ 
vation,  add  the  relevant  mean  hourly  slope,  add  a  correction  term  of  slope 
[Eq.  (8)]  that  is  a  function  of  deviation  of  foF2now  from  the  (4  day)  averaged  foF2. 
This  provides,  from  the  schemes  tested,  the  best  prediction  for  the  foF2  at  the 
next  hour. 

The  deployment  of  digital  ionosondes  with  automatic  data  processing  capabili- 

g 

ties  (Reinisch  et  al  )  will  make  real  time  half  hourly  data  available  to  the  EA  oper¬ 
ators.  Wo  plan  to  analyze  the  proposed  scheme  for  use  of  half  hourly  real  time 
data,  which  should  improve  the  predictions  even  further.  A  real  time  test  using 
Goose  Bay  and  Argentia  data  will  be  initiated  in  the  near  future.  A  further  im¬ 
provement  in  prediction  of  the  radar  frequencies  to  be  used  requires  the  incorpor¬ 
ation  of  the  change  in  layer  height  during  sunrise  into  the  prediction  scheme.  We 
are  currently  investigating  available  data  bases  as  to  their  usefulness  for  this 
analysis. 
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